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Abstract 

We study the electronic transport in a double quantum dot structure connected to leads in the Kondo regime for both series 
and parallel arrangements. By applying a finite- U slave boson technique in the mean field approximation we explore the effect 
of level degeneracy in the conductance through the system. Our results show that for the series connection, as the energy 
difference of the localized dot levels increases, the tunneling via the Kondo state is destroyed. For the parallel configuration, 
we find an interesting interplay of state symmetry and conductance. Our results are in good agrement with those obtained 
with other methods, and provide additional insights into the physics of the Kondo state in the double dot system. 
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Transport properties of small structures such as 
atoms, molecules and quantum dots have received a 
great deal of attention in recent years. Electron co- 
herence at these nanoscales bring out interference 
effects that are inherent to the system and are ex- 
plorable in experiments [1]. Confinement of electrons 
in these structures results in strong Coulomb effects. 
The combination of Coulomb interaction and strong 
coupling to reservoirs, for example, gives rise to inter- 
esting spin-singlet correlations observed at energies of 
the order of the Kondo temperature Tk, typically in 
the sub-kelvin regime. This Kondo effect has been ex- 
tensively studied theoretically and experimentally in 
different quantum dot systems [1,2,3]. Rich behavior is 
anticipated when coherence and correlations compete 
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in different geometries, and we address some of these 
issues in this work. 

Although theoretical studies have been presented 
utilizing different approaches [4,5,6,7], the physical 
richness of these systems deserves more attention. 
In this paper we focus our efforts on describing the 
transport properties of the double quantum dot sys- 
tem connected to leads either in series or in parallel 
while in the Kondo regime. We apply a finite- U slave 
boson mean field technique [8], which has been shown 
to capture the interesting physics in this regime. We 
find that the strong interdot coupling either directly 
or through the leads, results in interesting sub- and 
super-tunneling transport regimes, associated with the 
characteristic bonding and anti-bonding states of a 
diatomic molecule. It is fascinating that these features 
are projected into the Kondo regime, and result in 
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Fig. 1. Schematic of a double quantum dot coupled to leads. 

unusual interference effects that should be explorable 
in experiments with quantum dots. 

Figure 1 shows a schematic representation of our 
system. The transport properties will be studied by 
calculating the conductance between sites —1 and 1. 
We describe the system by a two-impurity Ander- 
son Hamiltonian H = H dot + Hi eads + H dot -i eads , 
where H dot describes the physics of the quantum dots 
(QDs) before they are coupled to the leads; Hi eads 
describes the two semi-infinite leads, and H dot -ieads 
establishes the couplings. We model the leads in 
the tight-binding approximation, which can be de- 
scribed by Hi eads = Y,a,i<-i (t c i* c i-if + H - c -) + 

J2a,i>i (^ c L c *+i^ + H.c.^j . In the finite-?/ slave boson 
approach, one enlarges the Hilbert space of the dots by 
introducing a set of boson operators, d (ej), pi a (pj^), 
and di (dj), which project onto the empty, singly and 
doubly occupied electron states respectively. Within 
this approach we can write the Hamiltonian H dot as: 

H dot = (£i + eV g )clc ia + U d Ui 

i=a,p i=a,/3 

er 

+ Y (tapzivCltvCpaZap + H.cj 

+ E (V^+E^qJ, (i) 

where 

zia = (l - d\di - p\ a V<^) ' (e\pi<r + P&d*j 

( + + \ -1/2 

X (l-ela-pl-p^ . (2) 
The constraints 



Pi = E p l p ^ + e * e * + d * di ~ 1 = ( 3 ) 

cr 

and 

Q ia - c\ a c i(7 - p\ a Pia - d\di - (4) 

are enforced in the problem through their respective 
Lagrange multipliers in Hdot, and A^, in order to 
eliminate unphysical states. Finally, 

H dot -i eads = + Y (tiJ z L c h c Jv + H - c -) ( 5 ) 

a,j = a,{3 
J = -l,l 

allows electrons to travel through the QD region. 

In the mean field approximation we replace all the 
boson operators by their expectation values. This pro- 
cedure results in a convenient non-interacting Hamil- 
tonian for the quasi-electrons with effective energy ei + 

(2) 

\ ia . The effective Hamiltonian is a function of 14 free 
parameters which are determined by the minimization 
of the free energy (H) with respect to all of them. The 
information about the strong correlation is however 
captured by these parameters. Utilizing the Hellman- 
Feynman theorem, the condition for minimal free en- 
ergy requires d{H)/dx = 0, where x runs over all the 
parameters. This results in a set of non-linear cou- 
pled equations that is solved numerically. The expecta- 
tion values appearing in these equations are calculated 
by Green's functions techniques. Expressions for the 
Green's functions can be obtained straightforwardly 
by the equation of motion method, using the effective 
mean- field Hamiltonian. From the Keldysh formalism 
the equilibrium conductance can be written as G = 
47T 2 t 4 pi(e F )pi(eF)\Gi 1 (e F )\ 2 , where pi(w) = pi (to) = 
p c (to) is the spectral density of states of the semi- 
infinite chain and G\i is the propagator that promotes 
one electron from site I to site 1 and cf is the Fermi 
energy. 

One can transform Fig. 1 into a series configuration 
if we take t a i = tp\ — 0. The lead bandwidth is D = At 
and the Fermi level is set to zero (ef = 0); hereafter we 
take the broadening T (= 7r£ /2 p c ) due to the coupling to 
the leads as the energy unit and show results for tem- 
perature T = 0. Figure 2 shows the conductance for 
the case of two dots coupled in series for several values 
of the energy difference Ae = ep — e a and for t a p = 1. 
For this relatively small value of t a p, the role of the 
direct interdot connection is only to provide a bridge 
between the two dots so that electrons can go through. 
Here, the level repulsion between the original dot levels 
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Fig. 2. (Color online) Conductance as function of gate volt- 
age V g for various values of Ae = ep — e a , for QDs in series. 
The energy unit is T; the parameters are U = 12.5 and 

*al = ^1 = * = °- 2t - 

is essentially negligible (within the level broadening), 
and the conductance has a maximum at the particle- 
hole symmetry condition (V g = —U/2). For the degen- 
erate case, Ae = 0, the dots are in the Kondo regime 
for the same V g , and the peaks in the density of states 
at the Fermi level provide a resonant tunneling condi- 
tion with unitary value. This limiting case agrees well 
with results obtained in previous works [5]. As Ae in- 
creases the dots reach the Kondo regime for slightly 
different V g values and the conductance is suppressed. 
Larger t a p values result in sizable level repulsion and 
two conductance peaks. 

Figure 3 shows the conductance as function of gate 
voltage for the parallel configuration for Ae = and 
different values of the interdot coupling t a p Notice that 
this configuration resembles the one used in Ref.[9] to 
study the persistent current of a quantum dot side- 
coupled to a ring. In this case however there are two 
quantum dots embedded in each arm of the ring. For 
the data displayed in Fig. 3, we have chosen U = 12.5 
and tpi — tpi — toil — t a i = t' — 0.25£. For t a p = 
we see the typical unitary limit conductance region 
around V g — —U/2. In this case we have both dots 
in the Kondo regime. The half-filled case results in a 
total spin S — 1, provided by the two electrons in the 
dots. One can distinguish two well-defined regimes on 
the conductance curves. The original flat region around 
V g = —U/2 for t a p = is split into two regions for 
larger interdot coupling. Let us analyze this behavior 
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Fig. 3. (Color online) Conductance as function of gate volt- 
age for various values of interdot tunneling t a p. The en- 
ergy unit is T; the parameters are U = 12.5, Ae = and 
t' = 0.25t. 

in more detail. 

As tap increases, level mixing and repulsion takes 
place, creating bonding and antibonding states in 
the dot system. These coherently-mixed states cre- 
ated between the two QDs have typical symmet- 
ric/antisymmetric configurations, resulting also in 
very different effective coupling to the leads. The an- 
tibonding state, for example, contributes a delta-like 
peak to the density of states (see Fig. 4) and is com- 
pletely uncoupled to the leads, due to destructive in- 
terference between the two paths. This interference is 
akin to the "subradiant" state discussed in the context 
of coupled radiant units, while here is more natural to 
call it the "sub-tunneling" state [10]. The contribution 
to the conductance, however, comes from the bonding 
state. The bonding/symmetric state has a much larger 
effective coupling to the leads, as in the "superradiant 
state", and is therefore termed the "super-tunneling 
state" [10]. The super-tunneling state results in the 
broad density of states feature seen in Fig. 4 at higher 
energies. 

It is the competition between the bonding and anti- 
bonding states that results in the conductance curves 
in Fig. 3. In the flat region present for t a p ^ 0, the an- 
tibonding state is found pinned around the Fermi level 
for a large range of gate voltage, as we can see in Fig. 4. 
In that regime, the contribution to the conductance 
comes from the tail of the bonding state that is slightly 
above the Fermi level. Since the bonding state is not in 
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Fig. 4. Density of states (left axis) as function of energy for 
a parallel configuration. Right axis indicate V g /U values 
for different curves (offset vertically). Same parameters as 
in Fig. 3 with t a(3 = 2T. 

resonance with the Fermi level, the conductance does 
not reach the unitary limit, until the antibonding state 
is filled and moves well below the Fermi level. Notice 
that the energy difference between bonding and anti- 
bonding states is determined by t a/ g. As the tail of the 
bonding state is increasingly occupied, the sharp anti- 
bonding state reaches the Fermi level. It remains there 
until it fills, while contributing nothing to the conduc- 
tance. Notice further that a change in the sign of t a (3 
produces a mirror reversal (left to right) of Figs. 3 and 
4. This is in agreement with the single-particle results 
in Ref. [11]. 

Breaking the degeneracy between dots, Ae / 0, sig- 
nificantly changes the conductance of the system. The 
conductance for the non-degenerate case is shown in 
Fig. 5. We have chosen here t' = O.lt and Ae = 0.075r, 
so as to have significant separation of the levels. We 
obtain a structure of three clear peaks, analogous to 
those obtained in previous works [6,7], where exact 
diagonalization plus embedding methods are used to 
study the Kondo regime. The three peaks can be un- 
derstood as follows. The peak around V g = is due 
to electrons travelling through the dot a with energy 
e a = V g — Ae/2, which is in the Kondo regime (S = 
1/2). As the levels shift with gate, the dot f3 with en- 
ergy e/3 = V g + Ae/2 becomes then accessible and the 
conductance starts dropping due to the destructive in- 
terference of electrons travelling through the two dif- 
ferent paths, totally vanishing the conductance at V g ~ 
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Fig. 5. Conductance for a parallel configuration for the 
non-degenerate case. The parameters are U = 12. 5r, 
Ae = 0.075, t' = 0.1 and t af3 = 0. 

— 0.217. The peak at the particle-hole symmetry posi- 
tion (Vg/U — —0.5) reaches unity when both dots are 
in the Kondo regime (S — 1). The following valley in 
the conductance has a similar interpretation as the first 
one. Finally, in the last peak (at V g ~ — 0.8U) the e a is 
well below the Fermi level and the contribution comes 
only from the dot /3 in the Kondo regime. 
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